Causality in Propagation of a Pulse in a Nonlinear Dispersive Medium 
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We investigate the causal propagation of the pulse through dispersive media by very precise 
numerical solution of the coupled Maxwell-Bloch equations without any approximations about the 
strength of the input field. We study full nonlinear behavior of the pulse propagation through 
solid state media like ruby and alexandrite. We have demonstrated that the information carried 
by the discontinuity, i.e, front of the pulse, moves inside the media with velocity c even though 
the peak of the pulse can travel either with sub-luminal or with super-luminal velocity. We extend 
the argument of Levi-Civita to prove that the discontinuity would travel with velocity c even in a 
nonlinear medium. 

PACS numbers: 42.65.-k, 42.50.Gy 



The propagation of a pulse of electromagnetic radia- 
tion through a linear medium depends critically on the 
dispersive properties of the medium. Sommerfeld and 
Brillouin investigated this problem in great detail 0. 
They discussed how the group velocity could be very 
different depending on whether one is working in the re- 
gion of anomalous dispersion or normal dispersion. Some 
of these results were verified 2J. Sommerfeld and Bril- 
louin also answered an important question - how infor- 
mation travels through the medium, fntcrest in this sub- 
ject has been revived [1 El since by using external 
laser fields one can produce a desired dispersion 0, @ • 
Such techniques have, in fact, led to the realization of 
ultraslow light. Based on an earlier suggestion of Chiao 
0, Wang et al. demonstrated the production of super- 
luminal propagation • The propagation of a pulse in a 
linear medium can also be understood from the interfer- 
ence of various Fourier components present in the input 
pulse JjJ. Stenner et al. were able to demonstrate the 
causal behavior in the propagation of pulses by studying 
the propagation of discontinuities in pulse shapes }12j . 

An important question is - how the discontinuities 
travel in a nonlinear medium? Needless to say that 
pulses in a nonlinear medium have been studied quite 
extensively since the early work of Mc-Call and Hahn 
[l3L H3. Ha. Hif . In this letter, we examine propaga- 
tion of discontinuities in a nonlinear medium. We choose 
solid state materials like ruby and alexandrite. Bigelow 
et al. studied pulse propagation in such materials and 
showed the possibility of both sub-luminal and super- 
luminal propagation 0, 0] . These authors specifically 
concentrated on the nonlinear regime. We use our ear- 
lier theoretical models [19( to study the propagation of 
the front of a pulse and its discontinuities in materials 
like ruby and alexandrite. We show how the argument 
of Levi-Civita [Ref.l,p.38] can be extended to prove that 
the front would travel with velocity c even in a nonlin- 
ear medium. This is true despite the fact that the group 
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FIG. 1: Schematic diagram of three level system for the ruby 
crystal. 



velocity could be either sub-luminal or super-luminal. 

Propagation of pulses with discontinuties in ruby — In 
order to develop an understanding of causality in pulse 
propagation through ruby, we consider the model as 
shown in Fig (Q. The figure shows the energy level of 
chromium ion's (Cr +3 ) in corundum (AI2O3) wherein we 
denote ground state level 4 A2 by \g), the 4 F2 absorption 
band by |ei), and the doubly degenerate levels 2A, E by 
\e2). An intense field E(z, t) — £(z,t)exp(ikz~iu!t) + c.c. 
couples states \g) and |ei) with Rabi frequency Q(z,t) — 
2d ls ■ £{z,t)/h, where d lg is the dipole matrix element 
and £ is the slowly varying envelope of the intense elec- 
tric field. Here we assume that the carrier frequency of 
intense field, lu, is in resonance with the frequency of the 
\g) < — > |e x ) transition. The density matrix equations for 
the system, as shown in Fig. are given by 0] 
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where p = (ejp|e ); i,j = 1,2. Ruby has very strong 
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relaxation effects where the dephasing rate I\ is very 
large compared to population relaxation rate r 2 (T\ 3> 
T 2 ). Typically T x ~ .f 5 GHz and T 2 - 35 Hz. Under the 
approximations Ti 3> T 2 , f2, we obtain p lg ~ and then 
the evolution equation for the ground state population 
can be written in the form 

n 2 

Pgg = i l - Pj-^rfrP gs - ( 2 ) 

In writing (0) we have used the moving coordinates 
i.e., the dot denotes the derivative with respect to 
2(t — z/c)r 2 . Note that under Ti ^> T 2 ,f2 we can set 
p lx « 0. In the slowly varying envelope approximation, 
we obtain the evolution of the Rabi frequency of an in- 
tense field as given by 

^ = --^&Pgg, & = SV^sat, ( 3 ) 

where the coupling constant ao = 4J\firu)\d lg \ 2 /cHTi and 
Q S at = v / 4rir 2 . In order to study the causality in the 
propagation of a pulse through ruby, we consider the 
shape of the pulse as incident on the input face of the 
medium given by 

Cl in = Cl°e- [T/a]2 t<T 

= t > T. (4) 

A sharply vanishing amplitude gives rise to the disconti- 
nuity in the pulse. Here 0° is a real constant indicating 
the peak amplitude of the pulse and a is the temporal 
width of the input pulse. Once the shape of the input 
pulse is specified, the working equations J2Jl and can 
be integrated numerically to investigate the evolution of 
the pulse. We consider all atoms to be initially in the 
state \g) at r = 0. The length of the medium L and 
coupling constant ao are chosen as 9 cm and 1.17 cm -1 
respectively. Fig. @ depicts the propagation of the in- 
put pulse through ruby crystal which is known to exhibit 
normal dispersion. To prove that the discontinuity of the 
pulse moves inside ruby with a velocity c, numerical in- 
tegration with high accuracy needs to be performed and 
the step size of the numerical integration to be taken of 
the order of 10 -10 in the time domain. This is due to 
fact that the discontinuity of the pulse takes only 0.3 ns 
to traverse a path length of 9 cm in the medium. We 
found that the difference in time between the discontinu- 
ity of the pulse at output and input ends of the medium 
is 65 x 10~ 9 in units of l/2r 2 by seeing the numerical 
output as given by table (P), which is not well resolved 
in Fig. J5J. Therefore, the discontinuity of pulse moves 
inside the ruby crystal with velocity c, although, it is ev- 
ident from Fig. (J2J) that the peak of the pulse propagates 
with a subluminal velocity. This conforms with causal- 
ity in the propagation of information through a nonlin- 
ear dispersive medium. We emphasize that in numerical 
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FIG. 2: The solid curve depicts light pulse propagation 
through a free space with a velocity c. The long-dashed and 
dashed curves show the light pulse propagating through the 
medium at different intensities. The discontinuity of the pulse 
moves with a velocity c irrespective of medium or free space 
and obeys causality. The temporal width of the pulse is 20 
ms and l/2r 2 = 4.45ms. 



TABLE I: Numerical solution of the coupled Maxwell-Bloch 
equations for the input pulse with peak amplitude f2° = 1 
and discontinuity occurring at 6.5 in units of 1/2F2 are given 
in this table. The path length traversed by the pulse through 
vacuum as well as the medium is 9 cm. The difference in time 
for the discontinuity of the pulse to traverse the medium is 
65 x 10 -9 in units of 1/2F2. All intensities are normalized to 
the peak values. 



Retarded time 
2(t - z/c)T 2 


Output through 
vacuum 


Output through 
medium 


6.4999999350 


0.945844116 


0.988440243 


6.4999999675 


0.945843694 


0.988440217 


6.5000000000 


0.945843273 


0.988440191 


6.5000000325 


0.945842852 


0.988440165 


6.5000000650 


0.945842430 


0.988440139 


6.5000000975 


0.000000000 


0.000000000 



simulations we have not used any kind of linearization on 
the nonlinear equations (0 and J2J. Note further that 
we work at intensities which are comparable to satura- 
tion intensities. The information travels in all cases with 
the velocity c even though the peak position depends on 
the power of the pulse. 

Propagation of pulses with fronts in alexandrite — To 
prove the causality of pulse propagation in a nonlinear 
anomalous dispersive medium, we consider alexandrite 
crystal where reverse saturation absorption is dominant. 
Reverse saturation absorption produces a narrow anti- 
hole in the susceptibility of the probe in the presence 
of a pump held which can lead to super-luminal prop- 
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FIG. 3: Four level model for alexandrite crystal. 

agation To investigate the propagation of the dis- 

continuity of the pulse through alexandrite, we refer to 
the ground state 4 A 2 by \g), the absorption bands 4 T 2 
and 4 Ti by |ei) and the level 2 E by \e2) as shown in 
Fig The intense pulse is defined by the electric 

field, E(z,t) = £(z,t) e - l M- fcz ) + c .c, which couples 
|ei) < — ► \g), and also drives |ea) < — ► |ea)- The density 
matrix equations for the model as shown in the Fig. 
are given by 



P 33 = (5) 

p 32 = -rap aa +M(fi„-PJ/2 

P la = -riP lB +ifi(p„-P u )/2 

Here we consider that the intense pulse is in resonance 
with |ei) < — ► \g) and |es) < — > leg), respectively. We can 
show that p 32 and p. vanish when the approximations 
ri,r3 ^> r2,f2; are taken into account. Under the same 
approximations, we obtain the equation for the evolution 
of ground state population as 



(i-p..)-|ftV 



(6) 



Note that p & P 33 ~ can be easily proved under 
ri,r3 ^> r2,f2. Therefore, the evolution equation for 
the slowly varying Rabi frequency of the intense field is 
given by 



on 

dz 



2 9 



2 V 



P..) 



(7) 



where ao and So gives the saturation and reverse satu- 
ration, respectively. Shand et al. have shown that, for 
an excitation wavelength of 457 nm, the excited-state 
Cl e 2)) absorption cross-section [a 2 ~ 4.05 x 10~ 20 cm 2 ) 




FIG. 4: The solid curve shows light pulse propagating at 
speed c. The dashed and long-dashed curves depict the prop- 
agation of the same pulse through an alexandrite crystal of 
length 9 cm with time delays — 51.7/4S and — 134^is respec- 
tively, whereas the discontinuity of the pulse moves with a 
velocity c. The width a and 1/2F2 are 500/is and 250/iS re- 
spectively. 



exceeds that of the ground state (\q) ) absorption cross- 
section (01 ~ 0.9 x 10" 20 cm 2 ) [2(J|. Following these 
experimental data, we estimate (evo/evo) ~ 4. The nu- 
merical integration of coupled Maxwcll-Bloch equations 
(JBJ and J7J) gives the evolution of the input pulse as stated 
by eq. Here we follow the same numerical integra- 
tion procedure as stated in the case of ruby. We obtain 
the time for the discontinuity of the pulse to propagate 
through alexandrite as well as vacuum to be 12 x 10~ 7 in 
units of 1 /2T2 by checking the output of the numerical 
result (table (dJ). This implies that the discontinuity of 
pulse moves inside the media with a velocity of light in 
free space c. We also found that the group velocity of 
the input pulse is in the range of superluminal velocity 
as illustrated in Fig. Q. Therefore, we present a very 
precise numerical simulation that proves that the discon- 
tinuity in a pulse of a electromagnetic radiation cannot 
propagate faster than the velocity of light c even though 
medium exhibits nonlinear anomalous dispersion. Con- 
sequently, the information carried by the discontinuity 
cannot propagate with superluminal velocity. It should 
be borne in mind that the used intensities in the nu- 
merical simulation are comparable with the saturation 
intensities and the solution of the coupled equations © 
and (2J do not use any kind of perturbation. Hence, we 
have studied the full nonlinear causal behavior of pulses 
through alexandrite. 

Finally we show how the argument of Levi-Civita, as 
quoted in the book of Brioullion pj , can be extended to 
prove causality for propagation in a nonlinear medium. 
Let us consider a plane electromagnetic field say x polar- 



4 



TABLE II: This table gives the result of the numerical inte- 
grations. The input pulse has peak amplitude f2 — 1 with the 
discontinuity at 4.2 in units of 1/2IV It shows that the time 
taken by the discontinuity of the pulse to propagate through 
medium as well as vacuum is 12 x 1CT 7 in units of l/2r 2 . 



GSA thanks R. Boyd and M.S. Biegelow for prelimi- 
nary data on their experiments. 



Retarded time 
2(t - z/c)V 2 


Output through 
vacuum 


Output through 
medium 


4.2000000 


0.9999999 


0.92618654 


4.2000003 


0.9999999 


0.92618647 


4.2000006 


0.9999999 


0.92618639 


4.2000009 


0.9999999 


0.92618632 


4.2000012 


1.0000000 


0.92618624 


4.2000015 


0.0000000 


0.00000000 



ized and travelling in the z-direction so that the electric 
and magnetic vector have the form E = xE(z,t)\ B = 
yB(z, t). The Maxwell equations then give us 



dE 

dz 
dB 

dz 



1 dB 
c~dt 

1 dE Air dV 
c dt c dt 



(8) 



where V is the nonlinear polarization of the medium. 
The medium is isotropic and hence the polarization is in 
x direction and is a function of z, t. Let e and h be the 
discontinuities in dE/dz and dB/dz at boundary of the 
wavefront. Note that at the boundary of the wavefront E, 
B and dV /dt are continuous. If v represents the velocity 
with which the discontinuity moves in the direction of its 
normal then the identity compatibility relations imply 
that 



discontinuity in 



-ve, 



and 



dE 

■ dB 

discontinuity m — — = —vn, 
dt 



(9) 



(10) 



This is seen from the equation E(z + Az,i + 
At) - E(z,t) = At dE/dt + Az dE/dz. Using com- 
patibility relations (|9I1()|I and Maxwell eqs. JHJ we get 
v = c. 

In conclusion, we have established by both numerical 
simulations and analytical methods that the front travels 
with the velocity of light in free space even in a nonlinear 
medium. 
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